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The treatment of the running QCD coupling in evolution equations is discussed. It is shown that 
the use of the virtuality of ladder (vertical) partons as the scale for QCD coupling in every rung of 
ladder graphs is an approximation that holds for DIS at large x only. On the contrary, in the small 
x region the coupling depends on the virtuality of s -channel (horizontal) gluons. This observation 
leads to different results for the Regge-like processes and DIS structure functions at small x. 

O : 

I. INTRODUCTION 

Of 

' In perturbative QCD, calculations of physical quantities such as scattering amplitudes or cross sections beyond the 
, Born approximation basically requires one to account for running a s effects. In fixed-order calculations, a s becomes 
' running in a straightforward way, by direct calculation of Feynman graphs. However, in various approaches, scattering 

amplitudes or cross sections are often calculated to all orders of the perturbation theory. Generally this means that 
fSJ ' some equations are constructed and solved and when such equations include the running of a s , the argument of a s 

is usually fixed by certain prescriptions which are based on different approximations. For example, for hard QCD 

processes exploiting DGLAP evolution (l), one takes 
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Si a s = a s (ki) (1) 
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in every rung of quark/gluon ladder, with k± being the transverse momentum of the ladder parton. As a consequence, 
in the evolution equations 



^ a s =a s (Q 2 ), (2) 

Ok _ 

O ' with \J Q 2 being the upper limit for integration over transverse momenta k± of ladder partons. In particular, when 
the DGLAP is used for calculation of the structure functions of the Deep Inelastic Scattering (DIS), —Q 2 is the 
i-^H virtuality of the photon. The dependence of a s given by Eqs (0j2|) follows from the results of ref. |^| for the treatment 
of a s in hard QCD processes. On the other hand, some important high energy QCD processes are Regge-like, rather 
than hard. In particular, high energy forward and backward scattering and DIS at small x are typical Regge-like 
processes. The use of the hard -like prescription Q of the QCD coupling for these Regge processes should have been 
justified somehow. Indeed, much recent work uses the argument of a s given by Eqs. (|l|, ^ without any comments. For 
example, in r ef § Eq. (0) is used for incorporating running a s into the BFKL equation though it is well-known that 
the BFKL is constructed especially for the Regge kinematics. Also, works Q suggest a s — a s (Q 2 ) for the polarised 
DIS structure function gi at small x. On the other hand, instead of exploiting directly Eq. (Q) refs [|| prescribe a 
modified parametrisation of a s , namely 

a s = a s (k 2 Jp), (3) 

where (3 is the parameter of the standard Sudakov representation of the ladder momentum k Q 

fc M = o% + /3p p + k± M . (4) 

Contrary to refs. [Q - Q, we have suggested a quite different treatment of a s in the evolution equations for the 
non-singlet structure functions at small x. It is interesting that when the argument of a s is discussed in ref j8|, the 
first assumptions about the argument look similar to Eqs. ([!]) though final result agrees rather with Q than with 



1 momenta p and q must obey 2pq S> p 2 ~ q 2 ~ 0, otherwise they are arbitrary 
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Eq. (|lj). In the present work we extend the treatment of a s suggested in M for the non-singlet structure functions, 
which are based on a quark ladder, to the the more complicated case of gluon ladders, forming the base for the singlet 
structure functions of DIS. We show that the treatment of a s suggested in refs. for hard processes and often used 
subsequently for description of different high energy QCD processes is an approximation valid only for hard QCD 
processes. This treatment cannot be used for the Regge kinematics. We show that the a s -dependence given by 
Eq. ([j]) is valid only for hard QCD processes where values of all Mandelstam variables are of the same order. In 
particular, for DIS it means that a s — a s {Q 2 ) (or a s (k 2 )) in evolution equations only when x ~ 1. On the contrary, 
the treatment of a s of Ref. |?J is valid for the evolution of DIS structure functions at small x and for other high 
energy QCD processes in the Regge kinematics. The paper is organised as follows: in sect 2 we discuss the treatment 
of a s in elementary parton processes. In sect. 3 we extend it to partonic ladders in the Regge kinematics and compare 
our results with results those of ref. 0]. Finally, sect. 4 contains our concluding remarks. 



II. INCORPORATING RUNNING a s INTO THE BORN SCATTERING AMPLITUDES 

In this section we discuss the argument of a s in elementary parton processes. In order to avoid unnecessary 
complications we regard all these partons as slightly off-shell. We specify this fact explicitly in Eq. ([l4|). Let us 
consider first the scattering amplitude M qq for the process of annihilation of a quark- antiquark pair into a quark 
-antiquark pair of the same flavour: 

?(Pi) + 9(P2) -> q(p'i) + ?(Pa) (5) 

in the forward kinematical region 

8 = (Pi +P2? « ~U= (p' 2 -pi) 2 » ~t = -{Pi-Plf . (6) 

In the Born approximation M qq orn is given by contributions of three Feynman graphs. In each such graph the 
intermediate gluon has momentum k, so that either k 2 = s or k 2 = u or at last k 2 — t. Such intermediate gluons are 
often addressed as s, u or t -channel gluons respectively. 

According to the Optical theorem 

crtot ~ s _1 3 s Af| f=0 , (7) 



so only the graph with an intermediate s -channel gluon in M qq orn contributes to the total cross section of the process. 
As our final goal is cross sections, we consider only this graph and denote its contribution by M qq orn . Adding the 
part of radiative corrections to M qq orn which makes a s running M qq orn transforms into 

Mqq = -47ra s (s) s + it ' 

where a s depends on s. In Eq. (|8|), we have dropped the factor An and the colour structure. The Feynman graphs 
corresponding to M qq are depicted in Fig. 1. The well-known BLM || arguments allow one to consider only quark 
loop contributions to a s and it is easy to show that the argument of a s in Eq. (||) is indeed s. To this end let us note 
first that at high energies and for space-like q obeying 

A 2 = A 2 QCD «-q 2 (9) 

a s is given by 

as{q2) =bH^W)' b=(UN-2n f )/(12n). (10) 



As is well-known, in Eq. (|10|) the contribution proportional to the number of colours, N comes from graphs with 
gluon loops whereas the contribution proportional to nj comes from the QED -like graphs where the gluon propagator 
is saturated only by quark loops. Obviously, the main, logarithmic contribution of quark loops is proportional to ln(s). 
As eventually quark and gluon loop contributions combine into Eq. (|Io|), the argument of a s in Eq. (||) is indeed s. 

Now let us repeat the above reasoning, replacing the quarks by gluons, i.e. let us consider the forward scattering 
process 
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g{pi) + gfa) -> gip'x) + g(p' 2 ) 



(ii) 



in the forward region Eq. (g). We denote the Born amplitude for the forward scattering of gluons by M^° rn and, 
similarly to Eq. (@), denote the Born approximation for the part of the forward amplitude with s -channel intermediate 
gluon by M^° rn (see Fig. 2). As is well-known, the Slavnov- Taylor identities state that accounting for radiative 

corrections to M g B g orn lead to the same dependency, a s = a s (s). Therefore both for the quark and gluon scattering, 
adding radiative corrections to the part of the Born scattering amplitude with an s -channel gluon transforms it into 



AI = -47ra s (s) 



R(s) 



(12) 



where we have dropped the color factors. R includes vertices and spin dependence. For example, for the scattering 
of on-shell quarks, R(s) — R q {s) — 2s. For the gluon scattering, R = R g depends on polarisations of gluons. We 
consider only its part, M, with s -channel gluon, instead of the whole scattering amplitude M, disregarding parts of 
M with the t and u -channel gluons because only M has non-zero imaginary part with respect to s. Eq. ( [l2| ) than 
leads to 



3 S M = SM = -3a s (s) 



^ttR{s) 



Ra s (s)ir5(s)(AnR(s)) 



(13) 



where we have used 3i?(s) = 0. We remind the reader that throughout this paper we assume all partons to be slightly 
off-shell. For the on-shell quarks, the second term in the rhs of Eq. (13) is zero. It is easy to see that a s (s), given by 
Eq. (O) when its argument is space-like, can be rewritten as 



a s (s) 



ln(s/A 2 



17T 



&ln(-.s/A 2 ) 6[ln 2 (s/A 2 ) +7T 2 ] ' 



(14) 



when s is time-like. Indeed, any scattering amplitude M{s) has a non-zero imaginary part in s when s is positive. In 
particular when M depends on s through logarithms, it means that 



M = M(ln(-s)) 



(15) 



when s is positive. However, Eq. (15) does not fix the phase of — s — se ltt} '. In order to fix <f>, one can use the fact that 
M does not have an imaginary part when s is negative. Then, the analytical continuation of M from positive s to 
negative s is defined through the upper path, because the cut in the s -plane when s > 0. In doing so, s acquires the 
factor e I7r leading to (j) = —tt. Thus, instead of Eq. ( |l5| ) one has 



M = M(ln(s) - m) 
when s is positive. Combining now Eq. (|l|) and Eq. @, we obtain 

An 2 R{s) 



3M 



6[ln 2 (s/A 2 ) +7T 2 



ln(s/A 2 )<5(s) 



(16) 



(17) 



Taking 3 s Af corresponds to cutting all intermediate states in the s -channel and summing up such contributions. 
The first term in the squared brackets in Eq. (17) corresponds to the case when the s -channel intermediate state is 



an on-shell gluon whereas the second term corresponds to many on-shell partons in the intermediate state. Using the 
(5-function in Eq. (|l^) when integrating the first term over s leads to a result which is singular in the infrared region. 
The reason for that is obvious: one can apply Eq. ( |Io|) for a s (s) only when s > A 2 . In order to fix this we replace 
S(s) by S(s — /x 2 ) in Eq. (J12|), with the infrared cut-off /*> A, arriving therefore at the expression 



3?M 



4ir 2 R(s) 



&[ln>/A 2 ) + 7r 2 



ln(s/A 2 )<5(s-/i 2 ) 



(18) 



for s -channel imaginary part of M qq and M gg . 

Finally let us consider the quark -gluon scattering, with both gluons off-shell, k 2 = (k 1 ) 2 < 



q{p) + g{k) -» q(p') + g(k') 



(19) 
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in the forward region (||). In the Born approximation the only Feynman graph with non-zero imaginary part in s is 
depicted in Fig. 3. Incorporating the radiative corrections to this graph in the same way as was done for quark-quark 
and gluon-gluon scattering and using the BLM approach ||, we conclude that for qg- scattering 

a s =a s (k 2 ), (20) 

which agrees with the treatment of a s suggested in ref. ||. Obviously, Q s a s (k 2 ) = and therefore 

instead of Eq. (pi). 



III. QCD COUPLINGS IN LADDER GRAPHS 

Now let us allow all partons in Eqs. (j^),(^l|) and (|l^) to be off-shell so that their virtualities cannot be neglected. It 
is obvious that letting the on-shell partons which were considered in sect. 2 be off-shell does not change our conclusions 
concerning the argument of a s in these scattering processes. Our formulae Eqs. (|lg|) , (pl|) also remain valid for off-shell 
parton scattering, save a change of R . However, letting the partons be off-shell converts the amplitudes we have 
considered into the rungs of partonic ladder graphs. Therefore we conclude that incorporating the running QCD 
coupling into the gluon-gluon and quark-quark rungs fixes the time-like virtuality of the intermediate s -channel gluon 
(the horizontal gluon line in Figs. 1-4) as the scale of a s . In contrast to this, the scale of a s for the quark-gluon and 
gluon-quark rungs is the space-like virtuality of the vertical gluon lines. 

Therefore the integrand of the quark ladder depicted in Fig. 4 contains n + 1 QCD couplings in such a way: 

a s (( P i - fci) 2 K((fci - k 2 ) 2 )a s ((k 2 - k 3 ) 2 )....a s ((p 2 + k n f ) . (22) 

When quark lines in Fig. 4 are replaced by the gluon lines, the integrand of such purely gluon ladder (we keep the 
same notations for the gluon ladder momenta as for the quark ones) contains the same succession of QCD couplings 
with the same arguments as Eq. (|2^). However, some interesting objects (e.g. for the singlet DIS structure function 
gi) are made of mixing of the quark and gluon rungs as shown in Fig. 5, so that the integrand corresponding to Fig. 5 
is proportional to 

a.(fc?)a.((fci) 2 )a a ((*3 - k,) 2 ) . . . a s (k 2 n )a 2 {{K) 2 ) (23) 
instead of the DGLAP succession^ 

a s (k 2 )a s (k 2 )a s (k 2 )....a 2 s (k 2 n ) . (24) 

A further difference concerning treatment of a s arises when those ladders are used for calculation of cross sections 
(e.g. DIS structure functions), which contain the imaginary part, 3 S A of these ladders: as all s -channel gluons 
(the horizontal lines in Fig. 5) are time-like, the related a s have non-zero Q s a s too. Accounting for them makes the 
pattern much more complicated than the standard DGLAP approach. 

Ladder graphs incorporating all orders in a s are essential ingredient of the DIS structure functions. We demonstrate 
now how and when it is possible to justify the DGLAP-like argument (|IJ) of a s for them. To this aim let us suppose 
that some DIS structure function <E> obeys a Bethe-Salpeter equation. The part of this equation involving the s 
-channel cut of the ladder virtual gluon is represented by the graph in Fig. 6. For example, it corresponds to the 
DGLAP equation or (after replacing the photon lines in Fig. 6 by gluon ones) - to the first term of the BFKL equation 
with running a s . After simplification of the spin structure, its contribution, $(s,Q 2 ) can be written as follows : 

^ 2 f d*k <l>(( g + fc) 2 ; Q 2 ,fc 2 ) 4fcj ^ a s (( P -k) 2 
(2tt) 4 (fc 2 ) 2 (q + k) 2 ^ S * (p-k) 2 



$(.s, <?) = (p + qY ,„ , 4 ' ; ,' n ± Z s 4tt * y n ( 2 5) 



where we have used the standard notations: q stands for the virtual photon momentum, Q 2 = —q 2 > and k± is 
transverse to the plane formed by q and p. The intermediate partons with momenta k in Fig. 6 can be either quarks 



2 for the ladder parton momenta k\ 
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or gluons. The external parton with momentum p is assumed to be on-shell and we will neglect its virtuality, which is 
not important for our conclusions. As we are discussing the treatment of a Sl we have dropped unessential numerical 
factors from rhs of Eq. ( p5| ) and have used the same notation <!> both for the on-shell structure function in the lhs 
of Eq. (||^) and for the off-shell one in the rhs. $ is supposed to be a logarithmic function of s. In contrast to the 
conventional DGLAP or BFKL cases, the rhs of Eq. ( jga ) contains 3a s instead of a s . Also, the argument of a s is the 
virtuality of the horizontal gluon. However under certain conditions, one can rewrite Eq. ( p5|) in the form containing 
a s with the DGLAP argument k\. To show this let us use the Sudakov representation for k 

k^ = a(q + xp)^ + Pp^ + kj_^ . (26) 

where x = Q 2 /s, s = 2pq, and then introduce the virtuality m 2 = {p — k) 2 of the s -channel gluon as a new variable 
instead of a. As k 2 = -{(3m 2 + k\)/{\ - (3), (q + k) 2 = [{(3 - x)[s(l - 13) - m 2 } - k\(l - x)]/(l - (3), we arrive at 

1 f„„2^^2 k l ^(W-x)[ S (l-P)-m 2 ]-kl{l-x)] 



(^ - R\ 2 n,Jm 2 \ 

(27) 



(1-/3) 

, 2 ((3m 2 + kl)^ s(l~(3) 2 ^ a s (m 2 ) 



(1-/3) ) [(/3-x)[s(l -(3) -m 2 ] - k 2 ± (l -xj] m 2 ' 

The integration over (3 is supposed to yield logarithmic contributions. It is necessary for BFKL and important also 
for DGLAP when applied in the small- x region. To this accuracy the integration over (3 in Eq. ( |27j ) can be written 
as dp/ (3, the main contribution coming from the region 

1 > (3 > 0min , (28) 

with 

Pmin X -\- — 7z (29) 

s — m 



as the lowest limit of the integration. Therefore Eq. (|27|) becomes 

dm 2 $(/3(s - m 2 ) - kl(l - x), Q 2 , {13m 2 + k 2 ± )"j 



/3 



k l ^(^s{m 2 ) 



(s — m 2 ) [(3m 2 + k 2 ^ 2 V m 



When the dependence of (3 m i n and $ upon m 2 in the rhs of Eq. ( |30| ) can be neglected, one can interpret the 
integration over m 2 as a dispersion integral for a s : 

77 Jo m [(3m 2 + k 2 ] 2 ^\ m 2 / 

1 r dm\„ S ^M^P) « asi ~^ m (3D 
V mr I k. 



7T 



[(3m 2 + k 2 ± } 



leading to an explicit a s dependence. In particular this is possible when 4> in the integrand of Eq. ( J30| ) docs not 
depend on k 2 and when x is big enough to satisfy 

f3min ~ X > k 2 ± (l - X)/ S . (32) 

Therefore the argument m 2 of a s in this kinematical condition is 

- ki/x < m 2 = -k 2 x /(3 < -k\ (33) 



Then the integration over (3 is in the region i</3<1 and Eq. (|32j) is equivalent to the condition k\[\ — x) < sx = 
Q 2 which is exactly the DGLAP region. According to Eq. (|32|), as x must not be small, the dependence a s w a s {k 2 1 _) 
of Eq. (|3^) actually coincides with the parametrisation suggested in Ref || . 

On the other hand, when x — > 0, (3 m i n depends on m 2 : 

Pmin » (34) 
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and, as well known, $ in the rhs of Eq. (pOf) depends on k . This makes impossible the use of Eq. (31) for performing 



the integration over to 2 . However, to the leading logarithmic accuracy, when m 2 <C s and \k 2 \ 3> (3m 2 , one can 

still perform the integration in Eq. without using the dispersion relation, as: 

1 f°°, fcj ^/ as(m 2 ) \ _1_ p 

Wo [/5m 2 + fci] 2 ^V to 2 J~£rP' 



R = - — 3a s (m 2 ) = - / -j (35) 

7T / A 2 TO/ 6 J n 2^ + 7T^ 



where we have used Eq. (|l4j) for ?ya s . Let us note that, contrary to Eq. (|3l|), the argument of a s in Eq. j35| ) is 
time-like. Introducing the infrared cutoff p 2 we arrive at the result 



1 f 1 dB [ sf) dk 2 t , , \ , , 



/« (36) 



for x — > 0, which actually reproduces the first term of the BFKL for the case when the gluon with momentum 
p is (nearly) on-shell. When n 2 is dropped, the integration in Eq. ( |35| ) yields R = a s (k 2 _/(3) otherwise R = 
(l/7rfc)arctan[ln(fci//3A 2 )/7r] . 

Eq. ( |36| ) shows that one cannot drop the dependence on 8 in R. In other words, there is no factorization between 
the dependence on the logitudinal and transverse momenta at small x. 

A simple way to solve such equations (and Eq. (|2q ) in particular) is known from the Regge theory (see e.g. 
PI). The point is that the small- x region is actually the Regge kinematical region where one can use the Mellin 
transform to simplify equations similar to Eq. (p5|). However, first one should introduce the signature amplitudes 
= [M(s,Q 2 ) ± M{— s,Q 2 )}/2 so that the structure functions are proportional to (l/7r)9 s M^ ± ^. Then one can 
use the Sommerfeld- Watson transform in the asymptotical form: 



100 doj ( s 



lOO 



Af(±)(.s/TO 2 ) = / f- (JL) e (±) (o;)F( ± )(o,) , (37) 
2m V to / 



where = (e~ 47r "±l)/2 are the signature factors well-known from the Regge theory. Although the transform ( |37| ) 
looks similar to the Mellin transform, actually there is a certain difference between them. Indeed, contrary to the 
usual Mellin transform, the transform inverse to ( |37| ) involves SjM^': 

F ( - ± 1(uj) = — -/ d(s/TO 2 )(s/TO 2 )- 1 - w 3 ;j A/( ± )( S /TO 2 ) . (38) 



sin(7ro;) J 

This method can be easily used to solve the Bethe-Salpeter equation (p5|), with the argument of 5a s being time-like. 
An application of this method has been given in for calculating the nonsinglet structure functions at small x. The 
fact that one cannot use the parametrisation a s — a s (k 2 _) for the DIS structure functions at small x suggests in 
particular that in the generalisation of the BFKL equation, when accounting for running a s effects, the first term of 
the new kernel should contain a s with the time-like argument kj_//3 (at least in the region where kj_/{3 dominates 
over the virtualities of the external gluons) whereas in the second term (corresponding to the virtual contribution) a s 
has the space- like argument k\. We will consider the running a s effects for the BFKL in a forthcoming publication. 

IV. CONCLUSION 

We have shown that when a s is running in a general QCD process, providing a s with the argument k 2 ^ is not a 
general rule but an approximation. This approximation is valid for hard processes but fails for the Regge-like ones. 
It can be used in the DGLAP evolution equations when they are applied in the kinematical region of large x but 
cannot be used at small x. In this regime we have suggested, for both the quark-quark and the gluon-gluon rungs 
involving Qa s , that the argument should be the virtuality of the s -channel (horizontal) gluon. Being incorporated 
into evolution equations of the Bethe-Salpeter type, this leads to an effective coupling with the argument k\ / (3. When 



J A similar parameterization was used in Refs. [fj, [pi for the non-singlet structure function. 



G 



the virtuality of this gluon is time-like, 7r 2 -terms appear in expression for Qa s (see Eqs. (^),(|TJ)), due to the analytic 
properties of a s . These 7r 2 -terms might be quite important because of their big numerical value. In particular, they 
have an impact on the values of intercepts of the non-singlet DIS structure functions |t]] and are likely to be important 
for the singlet structure functions as well. However, when the starting point /xo of evolution is big enouph, ir 2 -terms 
can be easily neglected. It can be used for estimating fx$. In particular, we obtained @ /^o ~ 5GcV for the non-singlet 
structure functions. 
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Figures 




FIG. 1. The Born amplitude M qq , having imaginary part at s > 0, with running coupling included. 




FIG. 2. The Born amplitude M gg , having imaginary part at s > 0. 
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FIG. 3. The Born amplitude M qg , having imaginary part at s > 0. 
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FIG. 6. A graph for DGLAP evolution of a DIS structure function with account of running couplig. 
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